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We study the interplay between disorder and a quasi periodic coupling array in an external mag- 
netic field in a spin-| XX Z chain. A simple real space decimation argument is used to estimate 
the magnetization values where plateaux show up. The latter are in good agreement with exact 
diagonalization results on fairly long XX chains. Spontaneous susceptibility properties are also 
studied, finding a logarithmic behaviour similar to the homogeneously disordered case. 
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Since their discovery in 1984 the properties of 
quasi-crystals have been a source of sustained interest. 
Many theoretical efforts on Ising models in Penrose lat- 
tices [H and XY Fibonacci spin chains ||-||] have re- 
vealed interesting magnetic orderings associated to the 
quasi-periodicity of these structures. Such kind of spin 
arrays have been found in recently synthesized rare earth 
(7^) ZnMg-7^ quasi-crystals (see e.g. [^) whose IZ ele- 
ments have well localized 4/ magnetic moments. The 
study of quasiperiodic ID chains has recently received 
renewed attention [^] and interesting properties have 
been elucidated. In^^ a system of spinless fermions in 
a quasiperiodic lattice potential was studied within per- 
turbation theory, where it was shown that its behaviour 
is different from both the periodic and the disordered 
cases: While in the case of a periodic potential one may 
have a metal-insulator transition only if the potential is 
commensurate, in the disordered case, the potential is 
relevant irrespective of the position of the Fermi level. 
In the quasiperiodic case, two different situations arise, 
depending on whether the Fermi level coincides with one 
of the main frequencies of the Fourier spectrum of the 
quasiperiodic potential or not. In the first case, the situ- 
ation turns out to be similar to the periodic case while in 
the second, at a perturbative level, the metal-insulator 
transition point is strongly modified. These predictions 
have been also verified numerically in Motivated by 
these studies we have recently analyzed the effect of an 
external magnetic field in a quasiperiodic spin chain, and 
found that the magnetization curve has a very interesting 
nature that could be predicted using a decimation pro- 
cedure and Abelian bosonization For the Fibonacci 
case, in particular, one can reproduce the main plateaux 
within the bosonization approach by approximating the 
quasiperiodic modulation by considering a subset of the 
main Fourier frequencies. From the experimental point 
of view, interest of quasiperiodic systems arise from arti- 
ficially grown qua si-p eriodic heterostructures |]l3|, quan- 
tum dot crystals JlJ] and magnetic multilayers [|15[. 

In this short note we go one step further to analyze the 
interplay between a quasiperiodic array of couplings and 
disorder in a XX Z spin chain in the presence of an exter- 
nal magnetic field. Using a simple decimation procedure 



we predict the appearance of plateaux in the magneti- 
zation curve at values of M which depend on both the 
quasiperiodicity and the strength of the disorder. As 
in the case of p-merized chains the presence of binary 
disorder results in a shift on the plateaux positions as 
a function its strength, while for the Gaussian case the 
plateaux are wiped out. We have also studied the ef- 
fects of a binary disorder on the double frequency XX 
chain studied in , where it turns out that the plateaux 
structure of the pure case disappears. We have checked 
this behaviour by studying numerically fairly long XX 
systems. We also predict a logarithmic behaviour of the 
susceptibility at low fields by extending the arguments 
in |1^] and have also verified this behaviour numerically. 
Other authors have also studied random spin systems re- 
cently @- @]. 

Let us consider the antiferromagnetic system 



Hqp — J ^ (1 + Cn) ( 5'n+l + 
n 



(1) 



where S^^S^^S^ denote the standard spin-^ matri- 
ces, in a magnetic field h applied along the anisotropy 
direction ( |A| < 1). Here, the coupling modula- 
tion is introduced via the Cn parameters defined as 
= COS (27rcjjyn ) , so quasi-periodicity arises 

upon choosing an irrational subset of frequencies cOj^ with 
amplitudes 6^ . 

Furthermore, the couplings Jn are randomly dis- 
tributed. Specifically, we consider a binary distribution 
of strength p {p = corresponds to the pure quasiperi- 
odic case while p = 1 corresponds to the uniform chain), 

P( J,) = p S{J, -U)^{l-p) 5{J, - J(l + en)) , (2) 

with en defined as above, along with a Gaussian disorder 
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These distributions, taken with 



same mean and variance, are built to enforce quasiperi- 
odicity. Thus, on average is a measure of the couplings 
quasiperiodicity. In what follows we assume that U is the 
smallest coupling 
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We will follow the decimation procedure as described 
in [|l^ to obtain the value of the magnetization for the 
main plateaux. In our problem (which is at T = 0) the 
energy scale is provided by the magnetic field, and in 
order to compute the magnetization, decimation has to 
be stopped at an energy scale of the order of the mag- 
netic field. We assume that all spins coupled by bonds 
stronger than the magnetic field form singlets and do not 
contribute to the magnetization, whereas spins coupled 
by weaker bonds are completely polarized. The magne- 
tization is thus proportional to the fraction of remaining 
spins at the step where we stop decimation. This simple 
argument happens to apply well to the binary distribu- 
tion, provided the energy scales of the involved exchanges 
are well separated. 

In studying irrational frequencies or other quasiperi- 
odic modulations, it is natural to analyze the case of 
the Fibonacci chain, a coupling array J a = J(l + (5), 
Jb = J (1 — ^) generated by iterating the substitution 
rules B ^ Ad^ndA^ AB [|, §, §, 0; with the dis- 
tribution (I) P(JO = p d{J, -U)^{l-p) S{J, - Ja,b)- 

We evaluate by decimation the magnetization of the 
widest plateaux in the strong coupling limit {5 ±1). 
There are two different cases to consider, according to 
6^—1^ i.e., Jb ^ Ja, and the opposite situation for 
1. 

Starting from saturation, in the first case the magnetic 
field is lowered until it reaches the value he ~ Jb at which 
the type-B bonds experience a transition from the state 
of maximum polarization to the singlet state. The mag- 
netization at this plateau is then obtained by decimating 
the B bonds. This yields: 



(M) = 1 



Nb 
Nt 



2(l-p) 
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where Nt = Nu + Na + Nb denotes the total number of 
bonds, Na,b the number of A and B bonds respectively 
and 7^ = {Na + Nb)/Nb- For a large iteration num- 
ber of the rules referred to above {Nt ^ oo), Na/Nb 

approaches the golden mean 7 = (^ ^~^^ ^ - In the p = 

limit, we recover the results in |^ and a non- vanishing p 
results in a shift of the position of the plateau. 

In the second case, J a ^ Jb, we have to distinguish 
two different unit cells since type- A bonds can appear 
either in pairs (forming trimers) or isolated (forming 
dimers). It can be readily checked that when lowering 
the magnetic field from saturation the first spins to be 
decimated correspond to those forming trimers. We then 
have a plateau (the nearest to saturation) at: 



(M)-L = 1-2 



Naa 
Nt 



1 
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where Naa refers to A pairs. The second plateau is ob- 
tained after decimation the type- A bonds, and then we 
must consider all the sequences J a between the others 
bonds. That gives for the second plateau: 



(M)2 = (Ml) 



^2 1 , 2 ^ 1 
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Again, we recover our results in ||] for p = 0. 

With this simple technique, one can predict the pres- 
ence and position of the plateaux, provided that there 
is a finite difference between the highest values of the 
couplings in the inequivalent sites. 

Since the decimation procedure applies for generic 
XXZ chains ||, we conclude that the emergence 
of these strong coupling plateaux is a generic feature, 
at least with an antiferromagnetic anisotropy parameter 
< A < 1. 

To enable an independent check of these assertions, we 
turn to a numerical diagonalization of the Hamiltonian 
(|l|) contenting ourselves with the analysis of the par- 
ticular case A = 0. This allows us to explore rather 
long chains using a fair number of disorder realizations 
(whose magnetization properties on the other hand, are 
self- aver aging). In Fig (^ we show respectively the whole 
magnetization curves obtained for various disorder con- 
centrations p = 0,0.2,0.4,0.6,0.8 and 1 averaging on 
5 X 10^ samples of L = /(18) = 2584 sites under the ex- 
change disorder ( ||), with 5 = 0.95 and U = 0.2. It can 
be readily verified that a set of robust plateaux emerges 
quite precisely at the critical magnetizations given by Eq. 
(^ for the plateau closest to saturation and by Eq. (||) 
for the second one. 

In a previous work ||] , we have observed that the mag- 
netization curve for the Fibonacci chain, could be well 
approximated by considering a rather small subset of the 
main frequencies in its Fourier spectrum. Here we study 
a two- frequency case, for a;i = 5/8 and 6^2 = 7/8, in 
the presence of disorder, where it is observed that the 
plateaux are erased even by small disorder (see Fig. (||)). 
It is interesting to observe that in contrast to the Fi- 
bonacci situation studied above, in which the plateaux 
structure is robust and just shifts with the strength of 
the disorder, here the plateaux seem to smear out even 
for a small value of p. 

It is important to stress that the derivation of our 
results for the quantization conditions Eqs. (||)-(||) rely 
strongly on the discreteness of the probability distribu- 
tion and would not to be applicable to an arbitrary con- 
tinuous exchange disorder. In fact, for the Gaussian case 
referred to above it turns out that no traces of plateaux 
can be observed. This is corroborated in Fig. (||), where 
we see that the plateaux structure is smoothed when the 
standard deviation is increased. Here the sampling was 
increased up to 4 x 10^ realizations though the length 
of the chain was reduced to L = /(18), as the CPU 
time per spectrum grows as . However, preliminary 
computations using larger chains yielded no substantial 
differences. 
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For homogeneously disordered chains, one can use the 
decimation procedure of along with the universality 
of the fixed point, to show that either for discrete or con- 
tinuous distributions the low field magnetic susceptibility 
behaves according to 



Xz oc 



(6) 



Following a simple argument based on random walk 
motion used in [|l9| , it can be readily shown that for A = 
(or XX chains), these arguments can be extended to 
the case of a disordered Fibonacci chain. It is interesting 
to note that the effect of the disorder is crucial since it 
changes the power law behaviour of the free Fibonacci 
chain obtained in Q to a logarithmic one. This can be 
clearly observed in the insets of Figs. (||), (H) where the 
validity of these arguments seem to apply over more than 
two decades. 

To summarize, we have studied the effect of disorder on 
the plateaux structure in quasiperiodic XX Z chains un- 
der an external magnetic field. By means of a simple real 
space decimation procedure we found the values of the 
magnetization for which the main plateaux emerge, Eqs. 
(||)-(||). This was tested by numerical diagonalizations of 
large XX chains finding a remarkable agreement with the 
quantization conditions in a variety of scenarios. Since 
the decimation scheme applies for generic XX Z chains 
JTsf , we conclude that the appearance of these plateaux 
is a generic feature, at least with an antiferromagnetic 
anisotropy parameter < A < 1. This issue still awaits 
numerical confirmation on sufficiently long chains using 
state of the art methodologies such as density matrix 
renormalization group [Q. Finally, we have also stud- 
ied the low magnetic field susceptibility which exhibits a 
clear logarithmic behavior, Eq. (^). We trust this work 
will convey a motivation for both experimental and nu- 
merical studies. 
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FIG. 1. Magnetization curves of modulated XX Fi- 
bonacci spin chains, immersed in disordered binary back- 
grounds of strength p after averaging over 5 x 10^ sam- 
ples with /(18) = 2584 sites, 6 = 0.95 , [/ = 0.2 and 
p = 0, 0.2, 0.4, 0.6, 0.8, 1 in ascending order. The left and 
rightmost lines denote respectively the pure uniform and pure 
Fibonacci cases. 
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FIG. 3. Magnetization curves of modulated XX Fibonacci 
spin chains, immersed in Gaussian exchange distributions, af- 
ter averaging over 4 x 10^ samples with /(18) = 2584 sites, 
8 — 0.95 and increasing standard deviation from left to right 
(note that the leftmost is practically the pure Fibonacci case). 
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FIG. 2. Double frequency magnetization curves of the XX 
chain for a;i = 5/8 and uj^ — 7/8 with amplitudes = 0.2 
and (52 = 0.3, with [/ = 0.1 and 10^ spins over 100 sam- 
ples, immersed in disordered binary backgrounds of strength 
p = 0, 0.2, 0.4, 0.6, 0.8 and 1 in ascending order 



h/J 

FIG. 4. Magnetic susceptibility of modulated XX Fi- 
bonacci spin chains, immersed in a Gaussian exchange dis- 
tribution [g — 1), after averaging over 4 x 10^ samples with 
/(18) = 2584 sites, b — 0.95. The inset show the susceptibil- 
ity behavior at low magnetic fields which follows closely the 
logarithmic regime predicted in the text. 
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FIG. 5. Magnetic susceptibility of modulated XX Fi- 
bonacci spin chains, immersed in a binary exchange disor- 
der of strength p = 0.6 averaged over 5 x 10^ samples with 
/(18) = 2584 sites, 6 = 0.95. The inset show the suscepti- 
bility behavior at low magnetic fields that like the Gaussian 
disorder follows closely the logarithmic regime. 
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